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Abstract 



Functionally Graded Materials are inhomogeneous elastic bodies 
j^ , whose properties vary continuously with space. Hence consider a half- 

space (x2 > 0) occupied by a special Functionally Graded Material 
made of an hexagonal (6mm) piezoelectric crystal for which the elas- 
tic stiffness C44, the piezoelectric constant eis, the dielectric constant 



O ! €11, and the mass density, all vary proportionally to the same "in- 

^ ' homogeneity function" /(X2), say. Then consider the problem of a 

piezoacoustic shear-horizontal surface wave which leaves the interface 
(x2 = 0) free of mechanical tractions and vanishes as X2 goes to infin- 
ity (the Bleustein-Gulyaev wave). It turns out that for some choices 

^f ■ of the function /, this problem can be solved exactly for the usual 

boundary conditions, such as metalized surface or free surface. Sev- 
eral such functions f{x2) are derived here, such as exp(zb2/3x2) (/? is 
a constant) which is often encountered in geophysics, or other func- 

t^ . tions which are periodic or which vanish as X2 tends to infinity; one 

final example presents the advantage of describing a layered half-space 
which becomes asymptotically homogeneous away from the interface. 

^ . Special attention is given to the influence of the different inhomogene- 

H I ity functions upon the characteristics of the Bleustein-Gulyaev wave 

(speed, dispersion, attenuation factors, depth profiles, electromechan- 
ical coupling factor, etc.) 



1 Introduction 

The wireless communication industry (mobile phones, global positioning sys- 
tems, pagers, label identification tags, etc.) fuels most of the current mass 
production of Surface Acoustic Wave devices (more than 1 billion units/year) 
where SAW-based interdigital transducers are used as high-frequency filters. 
In the race for miniaturization, devices based on Bleustein-Gulyaev waves 
(pure shear-horizontal mode) technology have proved more apt for down- 
sizing than those based on Rayleigh waves (two- or three-partial modes) 
technology, according to Kadota et al. (2001). 

Can the so-called "Functionally Graded Materials" , whose properties vary 
continuously in space, be used to improve the efficiency of Bleustein-Gulyaev 
waves? For a 6mm piezoelectric homogeneous substrate, the classic solution of 
Bleustein (1968) and Gulyaev (1969) is quite simple to derive; for a. function- 
ally graded substrate, the corresponding wave solution is in general impossible 
to determine analytically. In order to make progress, and with a view to use 
the eventual results as benchmarks for more complicated simulations, this 
paper strikes a compromise between these two extreme situations, and aims 
at finding in a simple way certain types of functionally graded substrates for 
which analytical Bleustein-Gulyaev type of solutions are easily derived. Such 
a task can be achieved by making the assumption that for the functionally 
graded material, the elastic stiffness C44, the piezoelectric constant 615, the 
dielectric constant en, and the mass density p, all vary in the same propor- 
tion with a single space variable. This assumption is often encountered in 
the literature, see for example the recent articles (Jin et al., 2003; Kwon and 
Lee, 2003; Wang, 2003; Chen et al, 2004; Kwon, 2004; Ma et al, 2004; Chen 
and Liu, 2005a, 6; Guo et al., 2005a, fe; Ma et al., 2005a, 6; Pan and Han, 
2005; Sladek et al, 2005; Sun et al., 2005; Feng and Su, 2006). It is a strong 
assumption, which can be envisaged to hold for C44, 615, en in certain con- 
texts (pre-stressed laminae (Cohen and Wang, 1992), elastic bodies subjected 
to a thermal gradient (Saccomandi, 1999), continuously twisted structurally 
chiral media (Lakhtakia, 1994), etc.), but is unlikely to hold for p as well. 
However, this assumption proves crucial for the derivation of analytical re- 
sults in terms of "simple" functions such as the polynomial, sinusoidal, and 
hyperbolic functions. Previous studies have indeed shown that if p behaves 
differently from the other material quantities, then analytical solutions of the 
shear-horizontal wave problem involve special functions such as Bessel func- 
tions (Wilson, 1942; Bhattacharya, 1970; Maugin, 1983), Hankel functions 



(Deresiewicz, 1962), Whittaker functions (Deresiewicz, 1962; Bhattacharya, 
1970), hypergeometric functions (Bhattacharya, 1970; Viktorov, 1979; Mau- 
gin, 1983), etc. (see the review by Maugin (1983) for some pointers to the 
wide hterature on the subject.); otherwise, numerical and approximate meth- 
ods are necessary to solve the problem, such as those based on Laguerre series 
(Gubernatis and Maradudin, 1987), on a combination of Fast Fourier Trans- 
forms and modal analysis (Liu and Tani, 1994), on Legendre polynomials 
(Lefebvre et al., 2001), on the WKB approximation (Liu and Wang, 2005), 
etc. 

In short, some generality is lost by taking p to behave in the same manner 
as the other quantities, but some simplicity and insights are gained, because 
the resulting exact solutions may serve as benchmarks for more realistic sit- 
uations, where for instance a perturbation boundary element method can 
be used (Azis and Clements, 2001). The governing equations derived in the 
course of this paper can also be specialized to the consideration of anti-plane 
deformations in the context of piezo-elastostatic problems, where the density 
plays no role and its eventual spatial variations need not be specified. 

The paper begins the analysis in Section 2 with the derivation of the 
equations governing the propagation of a Shear-Horizontal wave in the type 
of Functionally Graded Material just discussed. In Section 3, a change of 
unknown functions leads to the decoupling of the four first-order govern- 
ing equations into two separate pairs of second-order differential equations. 
For certain choices of inhomogeneity, the differential equations have con- 
stant coefficients and the consequences of such choices on the propagation of 
Bleustein-Gulyaev waves are fully analyzed and are illustrated numerically 
by two examples: one where the inhomogeneity is a decreasing exponential 
function, the other where it is an inverse quadratic function. The last two 
sections show that other inhomogeneity functions leading to explicit results 
can be generated, not necessarily by seeking differential equations with con- 
stant coefficients. Section 4 focuses on an inhomogeneity function for which 
the material parameters vary smoothly from a value at the interface to an 
asymptotic value at infinite distance from the interface. Section 5 presents 
one method, presumably among many others, to generate an infinity of in- 
homogeneity functions leading to exact Bleustein-Gulyaev solutions. 



2 A certain type of functionally graded ma- 
terials 

The Bleustein-Gulyaev wave is a shear horizontal wave, traveUing over the 
surface of a semi-infinite piezoelectric solid for which the sagittal plane is 
normal to a binary axis of symmetry. Now consider a half-space X2 ^ 
(say), made of a piezoelectric crystal with 6mm symmetry (see e.g. Royer 
and Dieulesaint, 2000) and with continuously varying properties in the X2- 
direction. Specifically, the elastic stiffness C44, the piezoelectric constant 615, 
the dielectric constant en, and the mass density p, all vary in the same 
proportion with depth X2'- 

{C44(X2), ei5(x2), eii(x2), P(X2)} = {C44, 6^5, 6°^^, p°} f {X2), (2.1) 

where C44, ejg, e"^, p° are constants, and / is a yet unspecified function of X2, 
henceforward called the inhomogeneity function. Without loss of generality, 
/ is normalized as /(O) = 1. 

Now take two orthogonal directions Xi, X3 in the plane 0:2 = such that 
the symmetry axis is along x^ and consider the propagation of a Bleustein- 
Gulyaev wave, traveling with speed v and wave number k in the Xi-direction. 
The associated quantities of interest are: the mechanical displacement com- 
ponent M3, the electric potential 0, the mechanical traction components ais, 
o"23, and the electric displacement components Di, D2- They are taken in 
the form 



{M3,0,^i3,^i}(a:i,X2,t) = {?73(x2),v^(x2),it,-3(a;2),irfj(x2)}e^'=(^^-^*), (2.2) 

where ?73, </), tj3, dj {j = 1, 2) are unknown functions of X2 alone, to be deter- 
mined from the piezoacoustic equations and from the boundary conditions. 
In the present context, the classical equations of piezoacoustics written 
in the quasi-electrostatic approximation, 

daij/dxj = p d^Ui/df, dDj/dxj = 0, (2.3) 

decouple entirely the anti-plane stress and strain from their in-plane coun- 
terparts. The anti-plane equations can be written as a first-order differential 
system, 
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and A7'2, K are the following constant symmetric matrices, 
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3 Some simple inhomogeneity functions 

In this Section, attention is restricted to some inhomogeneity functions for 
which the piezoacoustic equations turn into linear ordinary differential equa- 
tions with constant coefficients. 



3.1 Further decoupling of the piezoacoustic equations 

With the new vector functions u and v, defined as 



U{X2) = a//(x2)'u(x2), ^(Xa) = v{x2)/^/Jix2), 

the system (12.41) becomes 
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Now, by differentiation and substitution, an entirely decoupled second-order 
system emerges: 
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(3.3) 



and two simple ways of finding exact solutions for shear-horizontal wave 
propagation appear naturally. 
• Either (i) solve 
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(3.4) 



where Cq is a constant. Then the solution u to the second-order equation 
(I3.3l) i with (now) constant coefficients is easily found. Finally, u follows 
from (I3.ip i and v from the inversion of (I2.4l) i. 
• Or (ii) solve 
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(3.5) 



where Cq is a constant. Then the solution v to the second-order equation 
(13.31) 9 with (now) constant coefficients is easily found. Finally, v follows 
from (I3.IP 9 and u from the inversion of (I2.4p 9. 

Of course the two possibilities (13. 4p and (13. 5p do not exhaust the classes 
of solutions. The last section of this article shows how infinitely more inho- 
mogeneous proffies can be generated. 

Clearly now, if p is solution to (13. 4p . then — p is solution to (13.51) : so 
that if / is solution to (13. 4p . then 1// is solution to (13. 5p . The resolution of 
these two equations is straightforward, and the results are collected in Table 
1: according as to whether cq is positive, negative, or equal to zero (second 
column), several functions p{x2) (third column) and /(X2) (fourth column) 
are found. The inhomogeneity profile PO is common to the resolution of 
(13:1) and ([33]); proffies P1-P5 result from ([33D and P6-P10 from ([33]). The 
quantities (3 (inverse of a length) and 6 (non-dimensional) are arbitrary, so 
that the inhomogeneity functions /(X2) in P4 and P5 are essentially the same 
functions, and so are the inhomogeneity functions in P9 and PIO. 

Some of these inhomogeneity functions are often encountered in the geo- 
physics literature, such as the exponential function PO or the quadratic func- 
tion PI. Dutta (1963) used P2 for Love waves; Erdogan and Ozturk (1992) 
and Hasanyan et al. (2003) derived P0-P5 in a different (purely elastic) con- 
text; P6-P10 appear to be new, presumably because the preferred second- 
order form of the equations of motion is usually (I3.3|) i rather than (I3.3|) 9 (see 
(Destrade, 2001) for a discussion on this latter point.) 

The functions found present the advantages of mathematical simplicity 
and familiarity. Each of them however presents the inconvenience of describ- 
ing a somewhat unrealistic inhomogeneity, because each either blows up or 
vanishes as 0:2 —* 00, or blows up or vanishes periodically. These problems 
can be overcome by considering that they occur sufficiently far away from the 
interface, and by focusing on the near-the-surface localization of the wave. 

3.2 Exact solution 

Here the emphasis is on the complete resolution for the Bleustein-Gulyaev 

wave in Case (i) (proffies P0-P5). In Case (ii), the resolution is very similar, 

and the corresponding results are summarized at the end of this subsection. 

First, solve the decoupled, second-order, linear, with constant coefficients. 



differential equation fl3.3l) i for ii: 

u" + {PN2K - col)u = 0, (3.6) 

with a solution in exponential evanescent form, 

u{x2) = e"'""^!/", ^{q) > 0, (3.7) 

where t/° is constant and q is an attenuation factor. Then t/° and q are 
solutions to 

[eN2K-{co-eq^)l]U'' = 0. (3.8) 

The associated determinantal equation is the propagation condition, here 

[e{q' - 1 - {v/zf^Y - CoWiq' - 1) - Cq] = 0, (3.9) 

where v^ is the speed of the bulk shear wave in the homogeneous (/ = 1) 
material, given by 

P°v°t' = cI, + ell/el,. (3.10) 

The attenuation factors gi, g2 (say) with positive real part are 



qi = ^1 + Co/F - (i;/i;°)2, q^ = ^1 + cq/F, (3.11) 

provided the speed belongs to the subsonic interval 

0<{v/v°Tf <l+co/k^. (3.12) 

The smallest of these two quantities (gi) is indicative of the penetration depth. 
Here, the inhomogeneity affects the penetration depth in the following man- 
ners: for the exponential profile PO and for the hyperbolic profiles (P2, P3, 
P7, P8), the wave is more localized than in the homogeneous case (/ = 1); for 
the trigonometric profiles (P4, P5, P9, PIO), the wave penetrates further into 
the substrate; for the polynomial profiles (PI, P6), the penetration depth is 
the same as in the homogeneous case. Note in passing that the inequality 
(I3.12P puts an upper bound on the possible values of P for the trigonometric 
profiles P4, P5, P9, PIO (where cq = — /3^), namely: j3'^ < k'^, which means 
that the wavelength of the wave must be smaller than the wavelength of those 
profiles. These remarks are however preliminary and concern the behavior of 
the functions u and v with depth. The behavior of the wave itself is dictated 
by the functions u and v, see 03. ip . In particular, inequality 03.12p ensures 



that ■it(oo) = but, because u = {l/^/f)^i, not necessarily that u{oo) = 0; 
this latter condition must be tested a posteriori against each different form 

Now the constant vectors U^, U^ (say) satisfying fl3.8p when q = qi, 
q2, respectively, are easily computed and the general solution to (13.61) is 
constructed as: u{x2) = '-f^e~'^'^'^^^tJ-^ +'j2^~^'^^^'^U'^ where 71, 72 are constant 
scalars. Explicitly, 
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Then u follows from fl3.ip i as: u = {l/y/J)u, and v follows from the sub- 
stitution of this latter equation into the inverse of fl2.4l) i. which is: v = 
—ifN2~^u'. In the end, it is found that at the interface. 
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(3.14) 



Now the usual boundary value problems of Bleustein-Gulyaev wave prop- 
agation can be solved. For the metalized boundary condition, y:>{0) = and 
^23(0) = 0. These conditions lead to a homogeneous system of two equations 
for the set of constants {71,72}- That set is non-trivial when the follow- 
ing dispersion equation for the metalized boundary condition is satisfied for 
V = Vm (say). 
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(3.15) 

Here, the positive quantity x^ is the (bulk) transverse- wave electromechanical 
coupling coefficient. Recall that the classic Bleustein-Gulyaev wave is non- 
dispersive for a homogeneous metalized half-space. Its speed v°-^ is given 
by 



i<JvTr = 1 - x^ 



(3.16) 



Hence the effect of an inhomogeneity of the form found in Table 1 is readily 
seen from the comparison of the last two equations. 

For the free (un-metalized) boundary condition, ^23(0) = and ^2(0) = 
iA;eo</5(0), where eo is the permittivity of vacuum (see for instance (Royer and 
Dieulesaint, 2000, p. 310).) These conditions lead again to a homogeneous 
system of two equations for the set of constants {71,72}- The dispersion 
equation for the free boundary condition, linking the wave speed Vf (say) to 
the wave number is now: 
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(3.17) 



Comparison of Vf and v^ shows that Vf > Vm, whatever the choice of / in 
Table 1; so, by (13. lip , the wave penetrates deeper into the substrate when its 
surface is not metalized, as is the case for a homogeneous substrate. Recall 
that for a homogeneous half-space, the classic Bleustein-Gulyaev wave is 
non-dispersive for "free" boundary conditions, and that it travels at speed 
v°r given by 

KK)' = l-x7(l+eVeo)^. (3.18) 



Note that both Vf and Vm are such that (I3.12p is verified. 

In Case (ii) (profiles PO and P6-P10), it is found that the attenuation 
factors are still given by (13. lip , and that at the interface 
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(3.19) 



For the metalized boundary condition, the dispersion equation is now: 
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and for the free (un-metalized) boundary condition^ the dispersion equation 
is now: 
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(3.21) 
These equations could be rationahzed but this process might introduce spu- 
rious speeds. It can be checked that they coincide respectively with (13.151) 
and (13.171) for the exponential inhomogeneity function PO, and with (13.161) 
and (I3.18P for the homogeneous substrate. 

3.3 Examples 

Consider that the substrate is made of a functionally graded material for 
which the material properties at the interface X2 = ^ are those of a PZT-4 
ceramic (Jaffe and Berlincourt, 1965): cl^ = 2.56 x 10^° N/m^, e^j'g = 12.7 
C/m^, e°^ = 650 x 10^^^ F/m, p° = 7500 kg/m'^. The permittivity of vacuum 
is taken as: eo = 8.854 x 10"^^ F/m. 

When the substrate is homogeneous, the Bleustein-Gulyaev wave travels 
with speeds: v'^ = 2256.85 m/s and v°r = 2592.65 m/s, for metalized and free 
boundary conditions, respectively. The surface electromechanical coupling 
coefficient K^ is given by (Royer and Dieulesaint, 2000, p. 296), 

Kl = -4 — -^ ^ — — ^, 3.22 

o „,o2 I eo o2 ,,o2 ' V / 

the latter approximation being justified in the PZT-4 case. Here, Kg ^ 0.242. 
In the first example, the inhomogeneity function is decreasing exponential: 
/(a^2) = exp(— 2/3x2), /3 > (profile PO of Table 1). Then the mechanical dis- 
placement Us{x2) varies as: {1/ ^/J{x2}) exp —kqiX2 = exp —k{qi — (3/k)x2, 
and it is found here that 

^ (3.23) 




. P P e^i 
k^ k eo 

for metalized and free boundary conditions, respectively. Both quantities are 
clearly positive and the decay is secured. The dispersion equations (I3.15J) 

10 



and (13.171) give the wave speed in terms of the dimensionless quantity jS/k. 
The range for this quantity is chosen so that the inhomogeneity function 
decreases with depth in a slower fashion than the mechanical displacement 
for the metalized boundary condition — for the free boundary condition, 
the wave speed is so close to the body wave speed 



'f 



0.9999998^7 



that the displacement hardly decays at all. In other words, [3/k satisfies: 
2/3//c < gi — (3/k, where the right hand-side is given by fl3.23P i. This is 
equivalent to: (3/k < y/^ /{2^/l + x^) = 0.2015. Figure la shows the varia- 
tions of Vm (lower curve) and Vf (upper curve) with (3/k from (homoge- 
neous PZT-4 substrate) to 0.2. In this range, the inhomogeneity function 
has no noticeable influence on the speed of the Bleustein-Gulyaev wave with 
free boundary conditions, whereas it slows down significantly the Bleustein- 
Gulyaev wave with metalized boundary conditions, resulting in an increasing 
electromechanical coupling coefficient K"^ (Figure lb), from 0.242 to 0.324. 
For the second example, the inhomogeneity function is inverse quadratic: 
f{x2) = l/(/^^2 + l)^j (3 > (profile P6 of Table 1). Then the mechanical dis- 
placement f/3(x2) varies as: {1/ ^/f{x2)) exp —kqiX2 = {(3x2 + 1) exp —kqiX2- 
Here the decay is secured when (I3.12p is satisfied, which is equivalent to: 
V < vti the same condition as in the homogeneous substrate. For this pro- 
file, Co = and /'(O) = —2/3, so that the dispersion equations (I3.20p and 
(I3.2ip are easily solved. The metalized boundary condition gives: 
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and the free boundary condition gives: 
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For the purpose of comparison with the first example. Figures 2a and 2b dis- 
play the variations of the wave speeds with the dimensionless quantity (3/k 
over the same range ^ (3/k ^ 0.2. They show that the infiuence of each 
inhomogeneity functions is very much the same: here, the electromechanical 
coupling coefficient increases (from 0.242) to 0.310 instead of 0.324 for the 
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decreasing exponential profile. The present inverse quadratic inhomogeneity 
is however more satisfying to consider from a "physical" point of view, be- 
cause it decreases slower with depth than an exponential inhomogeneity, and 
it can thus describe a situation where the wave is confined near the surface 
while the material parameters (12.11) vanish at a greater distance. 

The next Section presents a third example of inhomogeneity function, 
this time yielding a profile for which the material parameters neither vanish 
nor blow-up with distance from the interface. 



4 An asymptotically homogeneous half-space 

Consider the well-known solution [C/(x2), V{x2)]'^ (say) to the piezoacoustic 
equations (12. 4p in a homogeneous substrate: 
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These functions satisfy (12. 4p when / = 1 that is, 

U' = iN2V, V = ik'^KU. 



(4.1) 



(4.2) 



Now seek a solution [u, vY to the piezoacoustic equations (12. 4p for the 
functionally graded material in the form: 
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jV + nv, 



(4.3) 



where Uq, Ui, Vq, Vi are yet unknown scalar functions of X2. By differentiation 
and substitution, it is found that if they satisfy 



then (12. 4p is satisfied. The choice 



0, v[ = 0, (4.4) 



Uo = l/\/7, ^'0 = V^, 



Ui = — — tanh5, Vi 
k 



/3 



(4.5) 



ktanh.6' 

(where /?, 6 are constants) takes care of (I2.4p i a.^. Then (12.4^ 9 ^^t both reduce 
to 

7)' - (/9/tanh5) = -(/5tanh(5)/, (4.6) 
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a solution of which is 



/(^2) 



tanh^(/3x2 + <5) 



(4.7) 



tanh 6 

This inhomogeneity function has the sought-after property of never vanishing 
(if 6 > 0) nor blowing-up as X2 spans the whole half-space occupied by the 
substrate. It describes a material for which the parameters C44(x2), ei5(a;2), 
eii(2;2)5 p{x2) change smoothly from their initial values C44, el^, e^i, p° at 
the interface X2 = to a higher asymptotic value {cl^, ejg, ej^, p°}/tanh^5, 
where 5 > is an adjustable parameter. The parameter /3 can also be 
adjusted to describe a not only a slow but also a rapid variation with depth, 
which would confine the inhomogeneity to a thin layer near the surface. This 
latter opportunity was excluded with the profiles of Section 3, because of 
their blow-up or vanishing behaviors. 

Now the substitution of / into 04.51) and then into the solution 04.31) . leads 
to the following expressions for the fields at the interface, 
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The dispersion equation for the metalized boundary condition is a quadratic 
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(4.9) 

At j3 = 0, the function / in 04. 7p is that of a homogeneous substrate (/ = 1), 
and this equation gives: 77 = x^, which, once squared, is 03.16^ . 

The dispersion equation for the free boundary condition is also a quadratic, 
now in ?7 = ^Jl — (vf/v^)'^: 
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At /9 = 0, this equation gives: 77(1 + Cii/^o) = id ■, which, once squared, is 

dSHD. 

Each dispersion equation is a quadratic giving a priori two roots: one 
tends to ?7 = (and so w = vt) as /5 —> (homogeneous substrate) and can 
be ruled out. 

For the third example, consider an inhomogeneous substrate whose prop- 
erties increase continuously according to (14 .71) from those of a PZT-4 ce- 
ramic at the interface X2 = to asymptotic values which are 10% greater 
(tanh5 = l/-\/l.l ~ 0.953). Figure 3 shows the variations of the inhomo- 
geneity function with depth, for several values of the parameter /3//c: clearly 
for fijk > 2, the inhomogeneity is confined within a layer near the surface 
whose thickness is less than a wavelength. Figure 4a shows the variations 
of Vm (lower curve) and Vf (upper curve) with j3/k from (homogeneous 
PZT-4 substrate) to 2. In this range, the inhomogeneity function again has 
no noticeable influence on the speed of the Bleustein-Gulyaev wave with free 
boundary conditions. However here the speed of the Bleustein-Gulyaev wave 
with metalized boundary conditions is always greater than in the homog- 
neous substrate, resulting in a smaller electromechanical coupling coefficient 
K^s (Figure 4b). 

5 More inhomogeneity functions 

The previous Section presented a method to derive an inhomogeneity function 
for which exact Bleustein-Gulyaev solutions are possible, but which did not 
rely on finding governing equations with constant coefficients as in Section 
3. That method is due in essence to Varley and Seymour (1988) (see also 
(Erdogan and Ozturk, 1992)) and it can be generalized to yield an infinity 
of such inhomogeneity functions. 

First, seek a solution [u, v]'^ to the piezoacoustic equations (12.41) in the 
form: 

/ -I hp—n / J kP~^ 

n=0 n=0 

(thus (14. 3 p corresponds to p = 1), where «„, Vn {n = 0, . . . ,p) are unknown 
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functions of X2. Differentiate once, substitute into fl2.4p . and use fl4.2p to get 

(5.2) 



and 

i («„ - /.„) v-'-'i + f ^ (!^ + „„ - /„„) v"-"-'> + v;v = 

n=l ^ ' 

(5.3) 
Both differential equations are satisfied when the following set of equations 
is satisfied, 

Mo = 1/V7> ^0 = \/7, Up = 0, t/ = 0, 

U'n-l/k + Un = Vn/f, v'n-l/k + f n = fUn- (5.4) 

Varley and Seymour (1988) found an infinity of / such that this set can be 
completely solved. Because there is little value in reproducing their deriva- 
tions, the reader is referred to their article for explicit examples. It suffices to 
notice that in general the solutions are combinations of trigonometric and/or 
hyperbolic functions, and that a great number of arbitrary constants are at 
disposal for curve fitting. For instance at p = 2, Varley and Seymour present 
at least 16 possible forms for /, each involving 5 arbitrary constants. As an- 
other example, it can be checked directly here that the polynomial function 

fix2) = i/3x2 + 1)'^ (5.5) 

where p is an integer, is also suitable; then, 

Un = bn{Px2 + lpP+^\ v^ = a„(/3x2 + 1)^'", (n = 0, . . . ,p), (5.6) 

where a„ and 6.„ are determined by recurrence from a^ = Bq = 1 and 

{p-n+l){p + n)an-i, hn = — ■ — a„, (n = l,...,p). (5.7) 



ikn p + n 
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Table 1. Some simple inhom,ogeneity functions 
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PO 


/32 


±2/3 


exp (±2/3x2) 
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2/5/ (/5x2 + 1) 


(/3X2 ± 1)2 


P2 


/5' 


2/3tanh(/3x2 + <5) 


cosh2(/3x2 + 5) 
cosh 5 
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/5^ 
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Figure 1: Influence of a decreasing exponential inhomogeneity function on 
the wave speed (free and metalized boundary conditions) and on the elec- 
tromechanical coupling coefficient. 
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Figure 2: Influence of an inverse quadratic inhomogeneity function on the 
wave speed (free and metalized boundary conditions) and on the electrome- 
chanical coupling coefficient. 
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Figure 3: Variation of an asymptotically homogeneous profile with depth for 
four different values of the dispersion parameter (3/k. 
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Figure 4: Influence of an asymptotically homogeneous profile on the wave 
speed (free and metalized boundary conditions) and on the electromechanical 
coupling coefficient. 
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